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ABSTRACT 

Performance improvement of a carbon dioxide gas-dynamic 
laser, due to regeneration, was studied theoretically. The 
study determined to v;hat extent preheating the nitrogen, 
before mixing it with the carbon dixoide and water, would 
effect available laser output energy. Combustion chamiber 
pressure and nitrogen inlet temperature were increm.ented 
by 100 psia and 300 K respectively. Due to the fact that 
little dissociation occurred at temperatures studied, no 
pressure dependence was found for the throat or combustion 
chamber temperatures. A m.ethod of preheating or regenerating 
the nitrogen is to use a gas-to-gas heat exchanger at the 
exhaust of the diffuser. Small signal gain was improved 
from 0.35 to 0.78 m~^. Limiting power extraction was 
increased by 42%. Optimum nitrogen temperatures were 
determined to be 1350 K for small signal gain and 2400 K 
for optimum power extraction. For typical operating con- 
ditions the fractional pressure loss across the heat 
exchanger was found to be 0.1 atmospheres. 
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I. INTRODUCTION 



The carbon dioxide gas dynamic laser has been considered 
a prime candidate for many high energy laser programs. 

One problem that has always plagued the high energy laser 
program is that of excessive weight and volume. These 
presently are the factors that are impeding the development 
of a tactical Airborne laser weapon system. A heat exchanger 
can be placed at the exhaust of the diffuser to preheat the 
nitrogen gas. This system significantly reduces the v/eight 
and volum.e of a given gas dynamic system while maintaining 
a constant level of energy output. 

The combustion chamber properties and nozzle throat 
conditions were determined by using a computer program 
developed by Sanford Gordon and Bonnie McBride at the NASA 
Lewis Research Center in Cleveland, Ohio. The Program and 
a detailed description of its use are contained in reference 
[1] , Although the program was written v;ith a rocket engine 
analysis in mind, it is adequate for computing equilibrium 
and frozen composition conditions in combustion gases used 
in gas dynamic lasers. 

After the computer determined conditions up to and 
including the nozzle throat, hand calculations were made 
using methods specified in reference [2] . According to 
calculations, the maximum available energy can be increased 



9 



by a factor of 5 to 10. The small signal gain as a function 
of nitrogen temperature reached a maximum value of 0.773. 
This maximum indicates a best operating temperature. 

The diffuser is used to recover the pressure in the flow 
stream and must be matched to the heat exchanger. 

A gas to gas heat exchanger is used to insure that the 
proper temperature difference is attained from the cold 
nitrogen side to the hot nitrogen side (see Figure 1) . The 
pressure must not be allowed to decrease to atmospheric 
pressure at mean sea level. This can be done by use of 
procedures outlined in references [3], (4], and [5], 
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Figure 1 BASIC ARRANGEMENT FOR A REGENERATIVE GAS 
DYNAMIC LASER WITH MOPA CONFIGURATION. 
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II. THEORY 



A. BASIC THEORY OF A GAS DYNAMIC LASER AMPLIFIER 

Carbon monoxide, hydrogen and oxygen are fed into the 
combustion chamber and ignited. The amount of each reactant 
is determined by the following reaction; 



9 ON 2 + 9C0 + H 2 + 5 O 2 -»■ 9 CO 2 + H 2 O + 9 ON 2 

After the reactant products reach equilibrium, nitrogen is 
mixed with the carbon dioxide and water in the mixing chamber. 
This yields a molar mixture of 90% nitrogen, 1% water and 
9% carbon dioxide. This is the composition required for 
operation of a gas dynamic laser amplifier. It is assumed 
that the ratio C /C of the mixture is 1.4 while the 
vibrational mode of carbon dioxide is frozen out at T^ . 

The nozzle array is designed to attain Mach 1.0 at the 
throat (region 3 in Figure 1) and Mach 4.0 or greater inside 
the cavity (region 4 in Figure 1) . For purposes of this 
thesis, all calculations have been done assuming that a 
Mach number of 4.0 exists inside the amplifier cavity. The 
technique used to determine the area ratios is described in 
detail in reference (6]. Since these calculations are well 
documented, only some of the final results will be mentioned. 
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From the equilibrium combustion temperature, the throat and 
cavity temperatures can be determined and the small signal 
gain can be calculated. 

The small signal gain is applicable only inside the 
cavity (region 4 in Figure 1) . Assuming that the oscillator 
and cavity are matched at a wavelength of 10.6 micrometers, 
the Output beam intensity is given by Equation 1. 

1 = 1^ e“'^ (1) 

a is the gain coefficient and I is the overall beam pathlength 
inside the cavity. £ is found from Equation 2. 

£ = NL (2) 



N is the nun\ber of times that the beam traverses the width 

L of the cavity (see Figure 1) . 

A measure of the power that can be extracted, defined 

in Equation 5.23 of reference [2] as the limiting power 

extraction, is a I where I is defined by Equation 3. 

' o o o J -1 



I 

o 



k(03-ei)QviB 






(3) 



The maximum available energy in the excited state of the 
gases is described in reference [2] and can be calculated 
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with Equation 4. Table II shows the result of calculations 
at incremental temperatures. 



E 



1 h\) 




(4) 



max 



m is defined as the average mass of the gas molecules, m 
is calculated with Equation 5 . 



ip is the mole fraction of the subscripted quantity. 

As a basis for determining the change in limiting power 
output, small signal gain and maximum available energy, an 
assumption was made that all gases were at a room tem.perature 
of 300 K, prior to combustion. The nitrogen temperature was 
then incremented and new calculations made at each value of 
nitrogen temperature. These new calculations are compared 
to those at 300 K in tabular and graphical form. 

B. SMALL SIGNAL GAIN 

Reference [2] outlines procedures for determining the 
small signal gain coefficient. The overall gain 1/1-^ may 
be determined from Equation 1 after has been found from 
Equation 6. Pages 5-12 of reference [2] shov;s the relation 
between a and a . 



m 




(5) 



o 
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-j(j+i) e 



a = 



1 o (2J+1) 
4tt tv^ T^/2 0^ 



fn(OOl) -n(lOO) , 

“TtS? ' 



'P 



CO. 



N 



( 6 ) 



The population inversion is determined by the quantity in 
brackets in Equation 6 . Population inversion is related to 
cavity and throat temperature by Equation 7, It is assumed 
that the vibrational temperatures of nitrogen and carbon 
dioxide are frozen at . Appendix B is a derivation for 
in terms of temperature and number density. 



rn(OOl) -n(lOO) , _ .e 

>• ill Kt <■“ 



-es/Ts 



- e 









"VIB 



(7) 



Here, the temperature is frozen at the nozzle throat and 
the system is operating so that the threshold condition is 
exceeded. That is, 63 / 62 . * "^his condition insures 

that the small signal gain will be positive. Q^-ib related 
to throat and cavity temperatures by Equation 8 . 



“ 6 ^/ T . ^ “® 2'^'^4 -2 -1 

Qvib = ) 



(8) 



C. MAXIMUM AVAILABLE ENERGY 

The maximum, available energy is dependent upon the 
wavelength of the laser and the vibrational temperature of 
the transition level for stimulated emission, E can be 
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calculated using Equation 4 or Equation 9. Equation 9 is 
the same as Equation 4 except that y has been substituted 
for V, 



E 



max 



1 ^ 
m X 





(9) 



D. HEAT EXCHANGER 

The heat exchanger is the key to success of this system. 

It takes energy from the exhaust gases at the diffuser and 
transfers it to the nitrogen gas. This process allows the 
nitrogen to be heated to an optimum temperature. Equation 
13 describes the pressure drop across the exchanger. Since 
this system is to operate at mean sea level and higher, the 
minimum pressure on the downstream side (region 7 of Figure 
1) of the exchanger must be greater than 1 atmosphere. 

Assuming that the nitrogen enters the exchanger at 300 K 
and exits the exchanger at a temperature of T^^ (optimal) , 
the total heat transfered is defined in reference [3] as 
Equation 10. 

must equal the heat given up by the mixture flowing through 
the heat exchanger. This is described by Equation 11. 

% = Qn, = Crsrr, [T^(HOT) - T (COLD)] (H) 

N ■‘m m pm m m 
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Therefore, Equations 10 and 11 can be combined to produce 
Equation 12 which describes the temperature drop across the 
heat exchanger. 



T 



5 




(?i) (V) 



m 



m 



pm 



[T^{HOT - T^(COLD)] 



(12) 



Here, is defined as the upstream temperature, (hot), 
and is the downstream temperature, (cold) . 

A more detailed analysis is in Appendix A. Equation 13 
is used in reference [4] to describe Ap. 



2 

Ap = -^5-Ji [(K +l-a^)+2{— - l)+f^ — - (1-a'^-K )— ] 
^ c V, ' A V, e V, 

i c i 1 



(13) 
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III. CALCULATIONS 



A. THE GAS DYNAMIC LASER CAVITY 

Calculations for the combustion chamber and nozzle 
throat were computed by using the NASA Lewis Computer 
Program, reference [Ij . The remaining calculations were 
>1 done by hand using the isentropic tables in reference [6] . 

For conditions in region 4, shown in Figure 1, Appendix 
A of reference [6] was used. At M = 4.0 the following 
relations are found: 



- Y 

^4 Y-1 2 

~ = (1 + M^) = 0.6586 X 10 

1 ^ 



=i = (1 + ^ M^) = 0.2381 

j-1 z 



A 

A^ " 



throat 



10.72 



Conditions in region 4 (Figure 1) can be determined from 
the following relationships: 



P 4 = (^)Pi = (0.6586 X 10 "^) P^ 
“^4 

T 4 = (^)T^ = (0.2381) T^ 



(14) 

(15) 
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mmm 










(10.72) A. 



(16) 



Ay, = (jT 

4 A 



-)A 



throat 



throat 



throat 



One model of the starting process uses a normal shock 
which develops inside the cavity (region 4 of Figure 1) . 
Therefore, the diffuser must be able to swallow this shock 
before the system will be stable enough to operate effec- 
tively, From Appendix B of reference [6] at = 4.0, the 
Mach number of the flow behind the normal shock is 0.4350. 
Before the diffuser can swallov; the shock, the flow must be 
accelerated to Mach 1 at the diffuser throat. Therefore 
from Appendix A of reference [6J at M = .4350, the following 
ratio is found: 



A 



A^ 



A 



A 



4 



diffuser 



1.5 



(17) 



Diffuser conditions change after the shock has been swallowed. 
Therefore, the only noteworthy condition is the fact that 
the ratio of diffuser throat area to nozzle throat area can 
now be determined. 



^diffuser 10.72 

^throat ■ ^ 



(18) 



By the methods described in Section 3.4 of reference [6], 
the Mach number at the diffuser throat is now 3.55 and the 
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normal shock is positioned just past the diffuser throat 
near region 6 of Figure 1. 



B. SMALL SIGNAL GAIN COEFFICIENT 

The small signal gain coefficient must be determined 
before the output intensity can be calculated. Combining 
Equations 6, 7, and 8 yields Equation 19. Equation 19 is 

explicit in terms of knovm throat and cavity temperatures. 

5 . . . 

The constant 1.6 x 10 is the result of inserting 

\ = 10.6 X 10“^m, J = 20, = 2.88 K, v = 39.5 x lO'^'^N JY? 

and T = 5.38 sec. 



-239.4 

T, 



= [1.6x10^] [- 



(TJ 



772 



] [- 



-3380 -1995 

-e 



T, 



“=Twr 



(l-e 



' )-^l-e 



T 



^77W 



^(1-e 



T 



3 )-" 



(19) 



Table 1 gives the values of nitrogen, throat and cavity 
temperatures used in calculating at the specified nitrogen 
temperatures. Figure 2 is a graph of versus temperature 
of the nitrogen. It is shown that an operating nitrogen 
temperature of approximately 1300 K will maximize the small 
signal gain at 0.773. See Table II for values of at 
specified nitrogen temperatures. 
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C. MAXIMUM AVAILABLE ENERGY 

The maximum available energy can be easily calculated 
by using Equation 4 after m has been determined by using 
Equation 5, 



- _ (.09) (44) + (.01) (18) + 
m — — 



(.90) (28) 



6.02 X 10^^ 

3380 



= 4.87 X 10 



-23 gm 



molecule 



E = (3.83 X 10 ) (e 
max 



- 1 ) 



-1 J oules 
Kgm 



These results are listed in Table II and Figure (3) is a 
plot of versus nitrogen temperature. 

D. LIMITED POWER EXTRACTION 

The saturation beam intensity, as defined by Equation 3, 
is a valuable aid in approximating the limiting power extrac- 
tion. Equation 3 must be used in conjunction with the 
following equations: 



-6 /T e /T -1 

= e ^ ^[(e ^ ^-1) + 2 






(e 



®2/^4 




( 20 ) 



e 

3 




) 



-6 J(J+1) 

2 _£ 

(2J+1) g *^4 
o W Tv^ T^/2Q^ ® 



(21) 



( 22 ) 
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Combining Equations 3, 8 , 20, 21 and 22 yields Equation 23. 



I 

o 



- 0 . 



-e. 



- 0 . 



k(03-0^) 



( 1 -e 



^4 -1 *^4 -2 "^3 

-^(1-e '^(l-e 



-1 



1 ^o 2J+1 

^3 ^ ^ “ T./TO" 

C ^ J- 



-0^J(J+1)/T^ 



(23) 



The product of Equation 19 and Equation 23 is Equation 24 . 



- 0 . 



- 0 , 



a I 
o o 

~ 7 ~ 



®3"®1 e 
^^2 ^4 



- e 



e_ ( — Tt P .+T .P . ) 
3 £3 14 4 4 



(24) 



2 

a I /P. is in Mks units. t.P. and T.P. can be found in 
o o 4 14 4 4 

Figures 2.8 and 2.10 of reference [2]. Table III shows the 

calculated results from Equation 24 and Figure 4 is a plot 
2 

of a^I^/P^ versus nitrogen temperature. 

E. HEAT EXCHANGER 

Using the optimum temperature as shown in Figure 4 
the hot nitrogen temperature is specified as 2400 K. This 
gives a temperature difference in the nitrogen of 2100 K. 

An energy balance will produce an equation relating the 
temperature difference from the hot to the cold side of the 
heat exchanger. From equation 10 and C^^^ equal to 0.39, 
the rate of heat transfer is 
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• • • RTIJ 

■3 = % V <“«> = 

Reference [7] contains values for C . From Equation 12, 

P 

the temperature difference is 

C ^ 

- T,^ = (1803) ~ 1617 K 

pm 

A more detailed analysis of the heat exchanger is in 
Appendix [A] . 

F. APPROXIMATE WEIGHT AND VOLUME 

Assuming one-second bursts of energy in the output, the 
amount of fuel needed for any pre-determined number of 
bursts can be calculated easiest at the nozzle throat where 
gas properties are well defined. From reference [6] , the 
mass flow rate is 

ra = UAp 

where U is the mass velocity, A is the area of the nozzle 
and p is the density of the gases. 

The limiting power per unit mass is 

a I 

= W 
m 
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where m = (1 sec) m . If W increases by some factor A and 

m 

m remains constant then a^I must have increased by A also, 

o o 

It is desired to maintain W constant and adjust the mass 
necessary to compensate for this increase in output. 



W(1 + A) 



(1 ^ 0^0 
m 



a I 
■ o o 

m 

TT+A) 



a I 
o o 

m ' 



m' 



m 

TT+aT 



This means that the mass necessary to maintain a constant 
output can be decreased by a factor of . Calculations 

have shown that A = 0.42 and the fuel necessary is only 70% 
of the fuel used without a regenerator. Therefore a savings 
in fuel weight and volume of 30% has been found. 
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IV. RECOMMENDATIONS AND CONCLUSIONS 



It has been shown that pre-heating the nitrogen gas 
prior to mixing it with the combusted products of hydrogen , 
oxygen, and carbon monoxide increases the gain of a carbon 
dioxide gas dynamic laser from 0.35 m ^ to 0.77 m~^. The 
limiting power extracted reaches a maximum at a temperature 
of 2400 K. This produces an optimum operating temperature 
and the limiting power extracted can be increased by a factor 
of 1.42. This implies an increase of 42% in output power 
while maintaining a constant mass flow rate. This means that 
a constant level of output power can be maintained while 
reducing the amount of fuels needed by 30%. This decreases 
the amount of storage volume by 30%. 

It is recommended that a further study be conducted to 
determine the effects that a heat exchanger would have on 
other fuels to produce a 10.6 micrometer wavelength with a 
possible extension into other wavelengths. 
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TABLE I 



TEMPERATURES AND PRESSURES OF SYSTEM 



TjJK) 





Tj (K) 


^4 (K) 


P^=P2(PSIA) 


P4(PSIA) 










500 


3.293 










600 


3.951 


300 


1076 


928 


256 , 


700 


4.610 


800 












_.900 


5.927 










1000 


6 . 586 










500 


3.293 










600 


3.951 


600 


1303 


1130 


310 I 


700 


4.610 


800 


5.268 










900 


5.927 










1000 


6.586 










500 


3.293 










600 


3.951 


900 


1535 


1338 


365 


700 


4.610 


800 __ _ 


5.268 










900 . 


5.927 










1000 


6.586 










500 


3.293 










^ _6oa .. 


3.951 


1200 


1774 


1553 


422 


700 


4.610 


800 


5.268 










900 


5.927 










1000 


6.586 










500 


3.293 










. 600 .. 


. 3.591 


1250 


1815 


1590 


432 


700 


4.610 


800 


5.268 










900... . 


__5.9.2L 










1000 


6.586 










500 


3.293 










600 


3.951 


1300 


1855 


1623 


442 


700 


4.610 


800 


5.268 










900 


5.927 










1000 


6.586 










■ 500 


3.293 










600 


3.951 


1350 


1895 


1662 


451 


700 


4.610 


800 


5.268 










900 


5.927 ^ 










1000 


6.586 










500 


3.293 










600 


3.951 


1400 


1935 


1694 


461 


700 


4.610 


800 


5.268 










_ . «00 


5.927. .... 










[IZMD 


6.586 
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TABLE I (Continued) 



T„ (K) 


T2 (K) 


T3 (K) 


T4 (K) 


P^=P2(PSIA) 


P4(PSIA) 










500 


3.293 










600 


3.951 


1450 


1975 


1730 


470 


700 


4.610 


800 


5.268 










900 


5.927 










1000 


^6.586 










500 


3.293 










600 


^ 3.951 


1500 


2014 


1772 


480 


700 


4.610 


800 


5.268 










900 


5.927 










1000 


6.586 










500 


3.293 










600 


3.951 


1800 


2246 


1991 


535 


700 


4.610 


SUIT 


~5TZ6'B 










900 


5.927 










1000 


6.586 










500 


3.293 










600 


3.951 


2100 


2463 


2204 


586 


700 


4.610 


800 


5.268 










900 


5.927 










1000 


6.586 










500 


3.293 










600 


3.951 


2400 


2662 


2404 


634 


700 


4.610 


800 


5.268 










5UD * 


027' 










1000 


6.586 










500 


3.293 










600 1 


3.951 


2700 


2845 


2590 


677 


700 


4.610 


5GD ' 


5 , Zb 3 










f 900 


5.927 










1000 


6.586 










500 1 


3.293 










600 ' 


3.951 


3000 


3017 


2760 


718 


700 


4.610 


800 


5.268 










900 


5.927 










1000 


6.586 
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TABLE II 



COMPUTED VALUES FOR E„^^ & a 

max 0 



T„ (K) 


(KJ/KG) ^ (KJ/lbm) 


aQ(nr') 


300 


10.3 4.67 


0.36 


600 


20.25 9.12 


0.57 


900 


33.3 15.1 


0.71 


1200 


49 22.2 


0.772 


1250 


51.9 23.6 


0.773 


1300 


54.5 24.7 


0.772 


1350 


57.7 26.2 


0.771 


1400 


60.3 27.4 


0.771 


1450 


63.3 28.7 


0.769 


1500 


66.8 30.3 


0.768 


1800 


85.9 38.9 


0.728 


2100 


105 47.6 


0.67 


2400 


124 56.3 


0.61 


2700 


143 64.9 


0.546 


3000 


159 72.1 


0.495 
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TABLE III 

COMPUTED VALUES OF a„I (^)(— ^)^ 

° ° ^ N/m^ 



T„ (K) 


T3 (K) 


T4 (K) 


' 


300 


928 


256 


0.543 


600 


1130 


310 


0.65 


900 


1338 


365 


0.746 


1200 


1553 


422 


0,827 


1350 


1662 


451 


0.844 


1500 


1772 


480 


0.8834 


1800 


1991 


535 


0.918 


2100 


2204 


586 


0.952 


2400 


2404 


634 


0.963 


2700 


2590 


677 


0.964 


3000 


2760 


718 


0.960 
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Figure 2 SMALL SIGNAL GAIN 
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Figure 3 MAXIMUM AVAILABLE ENERGY 
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Figure 4 LIMITING POWER OUTPUT 
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Figure 5 MIXING CHAMBER EQUILIBRIUM TEMPERATURE 
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APPENDIX A 






THEORY OF GAS-TO-GAS HEAT EXCHANGER 

As Stated on page 21 of reference [4] , the terms 
involving exchanger input and output flow are negligible. 
An inspection of Equation 13, shows that the assumption 



* 5 - ^ ^ - 1 ) » ^ 

c 1 1 1 

2 

can be used to develop Equation 4 ' . K^+l-o is associated 

2 ^2 

with the upstream flow while (1-a -K ) — is associated with 

e 

the downstream flow. 



2 

G v^ v_ V 

AP = -2^ [2(^ - 1) + f 



(4 ') 



By definition: 



V 

'^m 2 



Therefore ; 



AO - A , A, , ,f A 

- -2g- ‘TT'I — ^ + <7 S- 

^ 1 c c 



- 2 )] 



Now, the heat transfer area is needed. The derivation 
is based on the counter flow heat exchanger, from reference 
[7] , shown below. 
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T^(HOT) 


© ' 












1 

4 


' 


k 1 

< — 




‘ ©. 


i 


^ Tj^(COLD) 



The temperature variation as a function of distance 
through the heat exchanger is shovvTi below. 




q = hAAT (III) 

m SGc 



AT^ = logarithmic mean temperature difference 



h = heat transfer coefficient. 
A = transfer area 



Before A can be determined, AT is needed. Therefore 

m 



d(AT) = dT + dT., 
m N 
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where : 



= 

m 


change in hot gas temperature 


II 


change in cold gas temperature 



Now: 



dq = 


-Cm dT = C„ itL, dT„ 
m m m N N N 



Therefore : 



d(AT) 


^ dq + dq 

-C„ m„ ITL^ 

mm N N 


Also : 





dq = hAT dA 



Therefore : 




d(AT) 

"W) 


= h( K-) dA 

C„ iTL, Cm 

N w mm 



Integration of this equation leads to 

(COLD) ^ ^ 

Inlij, T ,,J0T) ' liv, “cm* 

5 N N N mm 
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As shown previously 



hA = q/AT 
^ ra 



and 



[T^ (HOT) - (COLD)) 

q/VTi^ = = 



AT 



m 



AT, 



m 



Now: 



T„-T., (COLD) T.-(HOT) - T., (COLD) - T_ + T_ 
In [?R — m — m Tvnr— j ~ 



^T_-T„(H0T) 
5 N 



AT 



m 



and 




[T^(HOT) - T^(COLD) - T^ + T^] 

T_ - T. (COLd) 

In [ ~ T - 1 

5 N(HOT) 



Now: 



A = 



"'mC„(Tc:-T7) p. u. A C (T. - T^) 
mmo/ DDCinb / 



hAT 



m 



hAT 



ra 



And the ratio A/A^ can be found as 



[p 



rUcC 

5 5 ra 



(Tg-T^) ]/[hATj 
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Assuming a 10% pressure drop across the exchanger, 
equal to 1.1 atmospheres and = 4500 °R, the Reynolds 
number can be calculated. 



R 



e 




P 



5 



y 



Assuming a tube diameter of .1 ft, -- 988 (ft/sec), 

— no 1 n~^ SJjUG - . * . n*”"^ SLUG 

p = 3.02 X 10 j and y “ 1.09 x 10 



ft' 



ft-sec ' 



R^ = 2.73 X 10' 
e 



From Figure 6-9 of reference [3] 



= 5.75 X 10^ 



''I'®''® '^ref 



and T j. are from Table a-5 of reference 
ref 



[3] . 



k 



(0.0337) 



4500' _ BTU 

I5W " ‘^""hr rt"«R 



h = 336 BTU/ (hr ft^ °R) 




Per Up C 
^5 5 m 

hAT 

m 




- ^ 7 ) 



12 



38 



0.02 



2 

With these quantities and ^ equal to 1.36 and f = 

^ 2 4 

-s, (from Figure 6-3 of reference [4] and G = 1.365 x 10 

lb^/(lbm ft sec^) ; 



AP 



2 

G V, 



'W 



0.1 atmospheres 
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APPENDIX B 



DERIVATION OF V 

c 



From Equation 3.7 of reference [2], is 



I N.o^C. = '*’N2 '’NjS2 ^ 



where 



\|) = the mole fraction 

a is the optical broadening collisional cross 
section . 

is the relative mean speed. 

♦cOj = -05 1-N^ = -90 -fHjO ' -Ol 

* 

From pages 3-5 of reference [2] 



-14 2 

o__ = 1.3 X 10 cm'" 

'-^2 

a., = .87 X lO"'-'* cm^ 

“2 

1 / o 

Q = .38 X 10 cm 
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From page 3.4 of reference [2] 



where 



Now: 



C. = 

^ 1/ Ttp 



y = for CO 2 collisions 



y = 



«COp«x 



for H 2 O and N 2 collisions, 



'CO, 



'N. 



'H2O 





K T^' 


\j 7T 


CN 

0 


r 


R T^' 




‘'^2 


r 


R T^' 



= 21 



(T 



M 

4 sec 



= 44 . 2 \|^ 



M 

sec 



I TT y 



= 40 \ T 



M 



4 sec 



«2° 



Therefore : 



= N X 10"^'*[{21) (.09) (1.3) + (44.2) (.9) ( 



+ (40) (.01) (.38) ] 



= 3.95 X 10 “^^ N 



.87) 
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